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Abstract 

Quantum random walks are shown to have non-intuitive dynamics, which makes them an at- 
tractive area of study for devising quantum algorithms for well-known classical problems as well 
as those arising in the field of quantum computing. In this work we propose a novel scheme for 
the physical implementation of a discrete-time quantum random walk using laser excitations of the 
electronic states of an array of quantum dots. These dots represent the discrete nodes of the walk, 
while transitions between the energy levels inside each dot correspond to the required coin oper- 
ation and stimulated Raman adiabatic passage (STIRAP) processes are employed to induce the 
steps of the walk. The quantum dot design is tailored in such a way as to enable selective coupling 
of the energy levels. Our simulation results show a close agreement with the ideal quantum walk 
distribution as well as modest robustness towards noise disturbance. 
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I. INTRODUCTION 



Quantum random walks represent a generalized version of the well known classical random 
walk, which can be elegantly described using quantum information processing terminology 
Despite their apparent connection however, dynamics of quantum random walks are 
often non-intuitive and deviate significantly from those of their classical counterparts j^. 
Among the differences, the faster mixing and hitting times of quantum random walks are 
particularly noteworthy, making them an attractive area of study for devising^fficient quan- 
tum algorithms, including thosepertaining to connectivity and graph theory 2|, y, |j] , as well 
as quantum search algorithms 5, ^. 

There are two broad classes of quantum random walks, namely the discrete- and 
continuous-time quantum random walks, which have independently emerged out of the study 
of unrelated physical problems. Despite their fundamentally different quantum dynamics 
however, both families of walks share similar and characteristic propagation behavior 

aaa. 

Strauch's recent work ol is the latest in a line of theoretical efforts to establishing a formal 
connection between the discrete and continuous-time quantum random walks, in a manner 
similar to their classical counterparts. 

There have been several proposals for implementing quantum random walks using a 



variety of physical systems including Nuclear Magnetic Resonance 
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14| . ion traps |l5|. classical optics 



16 



17 



19 



20 



21 



ll|, cavity QED 
, quantum optics 
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3^ as 



24i\ . optical multiports (25|, l26|, 1271] , optical lattice and microtraps [28, 
well as quantum dots (33j. 

In this paper we introduce another proposal for implementing the discrete-time or coined 
quantum walk on a line using a series of stimulated Raman adiabatic passage (STIRAP) 



operations 
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36( 1 on a single electron trapped in an array of quantum dots. An impor- 



tant advantage of our proposal is that it relies on well established and generally accessible 
experimental techniques which result in both high fidelity operations as well as a relative 
ease of scalability. To the best of our knowledge, the proposal of Solenov and Fedichkin 



33( 1 is the only other implementation to date which employs quantum dots, but unlike our 



scheme, it pertains to a continuous-time quantum walk on a circle. 

In what follows we present a brief overview of the coined quantum random walk (Sec. 
nil) and describe our proposal for the design of the quantum dot array and the sequence 
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of required STIRAP operations (Sec. IIIip . We then present a numerical simulation of the 
system's evolution (Sec. II VI) . including the effect of imperfect STIRAP operations. In 
the appendix we also demonstrate an efficient numerical technique for the optimization of 
STIRAP pulse parameters. 



II. COINED QUANTUM RANDOM WALK 

A one-dimensional quantum random walk consists of a walker hopping between nodes 
or quantum states \i) {—N/2 — 1 < i < N/2) assembled in a line. In the coined quantum 
walk, each state \i) further consists of two sub-levels or coin states labeled as | and 
I Unlike the classical case, the quantum walker has a complex valued distribution ijj 

over all the states, which remains undetected throughout the walk. Each step of the walk 
involves a coin flip, defined as a simultaneous unitary rotation 

/ cos(6') sm(e)e'^^ \ 
ysm{9)e"^^ - cos(^)e*('^i+'^2) J 

on the coin states of all nodes, followed by a conditional translation which shifts the walker 
in states | t, i) and | J., i) to states | |, z + 1) and | l,i — l) respectively. Hence for a quantum 
walker in an initial state i^{0), its state after n steps of the walk is given by ip{n) = U"- ipi^), 
where 

u = til ^ii c (2) 

is the overall evolution operator for a single step. A final probability distribution is deter- 
mined by collapsing the walker's wavefunction ip^n) at the end of the evolution. 
In this paper, we implement a modified evolution operator 

U = flC (3) 

which is the same as Eq. [2] up to a translation and relabeling of states. In other words we 
can define a mapping M : U i — > U where M = TLi L first relabels all the nodes according 
to L : \i) I — > \2i), followed by a translation T_i of the entire wavefunction one node to the 
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left, that is 



M U 



f_i L f C 



fii fli fl, C 

Til ^-Ii C 
U. 



(4) 



III. PHYSICAL IMPLEMENTATION 



To implement the quantum walk we use an array of quantum dots, depicted in Fig. [T^, 
where all the odd barriers have been significantly lowered creating pairs of coupled dots. 
Figure [T]d illustrates a more detailed structure of the pair, labeled as the walk quantum 
dot QD,^^;;. and the auxiliary quantum dot QD^„^. The essential feature of this design 
is that for low energies, the energy eigenstates of QD^y^;^ and QD^^^ are, to a very good 
approximation, spatially separable and the electron wavefunction is localized within the 
dot. For energies above their joint potential barrier however the two dots share common 
electronic states. 

The nodes of the quantum walk are mapped to the successive QD^^/^;^ along the array of 
quantum dots. As depicted in Fig. [21 the first two energy levels of QD^y^;;. encode the coin 
states I I) and | T) of the walk, the fourth energy level of QD^^^ represents an auxiliary 
state I A), Ee represents an excited state |e) well above the joint barrier between the two 
dots, and other states remain unoccupied throughout the walk. The quantum walk itself is 
represented by the propagation of a single electron wavefunction through the array of dots 
using a series of specially optimized 2- and 3-photon A STIRAP operations. 

The 2-photon STIRAP is used to perform the translation operation T|i. Here two laser 
pulses, pump P and Stoke 5, with angular frequencies VL^ and VLa respectively, couple the 
dressed states | t) and 1^4) via the intermediate state |e). By tuning the laser parameters 
and applying the two pulses in the counter intuitive sequence, one can achieve coherent 
population transfer between states | t) and \A) with almost perfect fidelity and without 
leaving any appreciable papulation residual in state |e) (See Fig. [5]). 

Likewise a pair of 3-photon STIRAP is used to perform the coin operation C. First 
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three laser pulses, PI and P2 and 5", with angular frequencies and respectively, 

couple the dressed states | J.), | t) and \A) via the intermediate state |e). A second 3-photon 
STIRAP is then applied in the reverse order. This procedure is shown to be capable of 
performing arbitrary rotations on the superposition state a\ |) |), independently of the 
initial amplitudes a and (3 jsT] (See Fig. E]). Therefore with careful optimization of laser 
parameters we can engineer a variety of coin operators C. 

Figure H] illustrates how the quantum walk can be performed. First a pair of 3-photon 
STIRAP operations perform a coin rotation C simultaneously on all coin states. A 2-photon 
STIRAP will then transfer all the | states to their corresponding \A,i) state. We then 
adiabatically raise all the odd barriers and lower the even barriers, virtually reversing the 
paring of the quantum dots. In this new arrangement every QD^^^ previously associated 
with the ith node is now paired up with (z + l)th QD^^^^^. Using a second 2-photon STIRAP 
we can now transfer \A, i) to \ + 1) which completes the implementation of T^^ operator. 
The barriers are then returned to their original setting and the process repeated for additional 
steps. The final quantum walk distribution corresponds to the probability distribution for 
detecting the electron inside each QD^y^;;. in the array of dots. 

An important consideration in the design of the quantum dots is the ability to perform 
selective addressing of states which are being coupled via STIRAP and to avoid all unwanted 
secondary excitations. Taking, for example, the pulse which is intended to couple the 
energy levels and Eg, the quantum dot energy structure should disallow a secondary 
excited state, say Ee' = Ei + Q^h to exist, as it would lead to the unwanted excitation of 
the El level. Similarly, to avoid leaking the electron out of the dot via ladder excitations, 
the energy structure should prevent the coupling of E^ to an upper energy level E^ + flih. 
What is attractive about our proposal is that experimentally this can be achieved without 
resorting to complex profiles for the quantum dot potential. In fact using simple square wells 
with dimensions given in Fig. [T]d we were able to produce the necessary energy structure 
depicted in Fig. [3l Assuming absorption line widths a < 1 meV, all superfluous excitations 
will be far off resonance and will not have any appreciable magnitude. 
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IV. RESULTS 



In order to simulate the evolution of the quantum walk in the array of dots, we first 
tuned our laser parameters to correctly perform the desired C rotation and translation 
operations. For our 2-photon STIRAP operations we employ pump and Stokes pulses with 
Gaussian envelopes Ss{t) and Sp(t) and parameterize them using their peak interaction ener- 
gies Sp and Ss, standard deviations ap and as, phase angles ap and as, and the time interval 
AT between the peak interaction energies. The STIRAP process can now be modeled by 
the time-dependant hamiltonian 



( ^.(t)e^"» ^ 
\^ Ep{t)e-'^^ y 



(5) 
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in the Raman resonance limit. 



constructed using the Rotating Wave Approximation 
The pulse parameters need to be tuned such that the resulting 2-photon STIRAP operation 
coherently transfers an entire population from state | t,^ to state \A,i) via state \e,i). We 
achieve this by optimizing Ep and AT, using a technique detailed in the appendix, while 
other parameters are fixed to any desired values. Figure [5] shows the time evolution of the 
dressed states | ^,1) and \A,i) and |e,z) under the application of the optimized 2-photon 
STIRAP, where we have set Eg = 1.5 meV, (Jp = (Xs = 4.0 ps and ap = as = 0, and the 
optimum Ep = 1.50 meV and AT = 5.87 ps. In order to achieve the second transition from 
\A, i) to I I, i+ 1) (after raising and lowing the alternate potential barriers) we simply reverse 
the order in which the pump and Stoke pulses are applied. 

In the double 3-photon STIRAP process depicted in Fig. [6l the pulse parameters need 
to be tuned to perform a unitary operation C on the coin states. As before we achieve this 
by optimizing E^ and AT after fixing the other parameters to any desired values. Setting 
Ea = 1.0 meV, a = 4.0 ps, = Opi = Pp = 0, ap2 = Psi = Ps2 = vr, api = Psi = vi", 
C(p2 = l3s2 = and as = 0, with optimum parameters E^ = 1.34 meV and AT = 6.12 ps, 
we obtain symmetric coins by simply varying jSp. Figure [6] depict the time evolution of the 
dressed states | |), | t), \A) and |e) under the action of the coin operators (7(^4, V^iV^) and 
C {11/6,^/2,^/2) for i3p = Y2 and /5p = '^/s respectively. We also obtain asymmetric coins like 
C{^/4, V2, -V2) by setting a^i = fSsi = it, ap2 = /3s2 = V2, «s = and /3p = '^/2. 

Following the control pulse optimization, we obtain the full 3x3 and 4x4 evolution 
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matrices corresponding to the 2- and 3-photon STIRAP operations respectively. We then 
use these to simulate the evolution of a single electron under the repeated applications of 
the pulse sequence outlined in Fig. m In Fig. [7] we have plotted the electron wavefunction 
after 100 steps, using optimized pulses corresponding to the translation operator T^-^ as well 
as two different coin operators C{'^/a^ and C*('^/6, '^/2, ^/i). The results are in excellent 

agreement with their respective ideal theoretical distributions. We also investigated the 
effect of noise disturbance and experimental uncertainty on the resulting distribution and 
demonstrated a relatively robust response against imperfect pulse parameters. Figure [8] 
shows a reasonable degree of fidelity after the introduction of white noise in the energy 
peak, phase, timing and the standard deviation of the laser pulses. 



V. CONCLUSION 



We have proposed a physical implementation of a discrete-time quantum random walk 
using the action of 2- and 3-photon STIRAP operations on an array of quantum dots. We 
demonstrated that our scheme reproduces the characteristic quantum walk probability distri- 
bution which remains observable after the introduction of modest experimental uncertainty 
in the laser excitations. 

Like many other proposed schemes however, our implementation of the quantum walk is 
essentially a wave interference experiment and does not involve any quantum entanglements. 
Such implementations come with a cost as the number of resources grows, at best linearly 
with the number of nodes required for the walk. Furthermore it is generally expected 
;hat almost all potentially useful applications of quantum walks such as search algorithms 



5| or element distinctness 



39|, stem from higher dimensional walks on general graphs. 



Nevertheless implementing one dimensional quantum walks is significant for carrying out 
feasibility studies of assembling such physical systems. 
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APPENDIX: Control Pulse Optimization 



Considering a time- dependant hamiltonian H{t) for the 2-photon STIRAP, its action on 
a three-level system can be determined by solving the Schrodinger equation 



H 



^2(0) 

V^3(0) ) 



(A-l) 



We do this by approximating H{t) using a series of time independent Hi over suitably short 
time steps which allows us to write the solution as 



MO) 



(A-2) 



where the evolution operator 



^-iHtSt/h^-iHt-i&t/h _ _ _ ^-iH25t/h^-iHi5t/h 

^ un U12 Ui3 ^ 



U21 U22 U23 
U31 U32 U33 



(A-3) 
(A-4) 



By carefully optimizing the pulse parameters we can achieve 



T 





(A-5) 



which is the ideal swap operation between states |1) and |3) via the intermediate state |2). 
When state |3) is initially empty, this amounts to a translation operation which coherently 
transfers an amplitude from state |1) entirely to the empty state |3) without populating the 
intermediate state |2). 

We achieve the optimization by first fixing Ss and phase angles dp, ag, ap and a^, and 
then varying Sp and AT in order to minimize the cost function 



Kt 




(A-6) 
(A-7) 



The exponentials ^-"^^tM/f^^ which have to be re-evaluated for every 
are efficiently and accurately computed using a Chebyshev expansion 

M 

= y^an(a)0n(^), 

n=0 



parameter variation, 



40 



4l| 



(A-8) 



where A = —iHt6t/h, a„(a) = 2J„(a) except for ao(a) = Jo{a), Jn{c() are the Bessel 
functions of the ffist kind, 0„ are the Chebyshev polynomials, and A/" is the number of terms 
in the Chebyshev expansion. To ensure convergence, the exponent A needs to be normalized 

as 

~ 9 A 

A = — , (A-9) 



where /imm and ^max represent the minimum and maximum eigenvalues of A. Chebyshev 
polynomials 0„ are efficiently evaluated using the recurrence relation 



,{A) = 2A^n-M)+(l)n-2{A), 



(A-10) 



and 



</.o(^) = l, <PM)=A- (A-ll) 

In practice, iterations are continued until the norm of the matrix exponential converges to 
the required level of accuracy. 

The 3-photon STIRAP process is similarly represented by 



Uc{t) 



where the evolution operator 



Uc{t) 



( ^i(O) \ 

V^2(0) 
^3(0) 

V ^4(0) / 



(A-12) 



uii U12 ni3 ni4 

U21 U22 U23 U24 

U31 U32 U33 

\ M41 M42 M43 M44 / 

This time pulse parameters can be optimized in order to to achieve 



(A-13) 



Mil U12 
U2I U22 



(A-14) 
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where (7 is a desired unitary coin matrix given by Eq. [H 

As before, this is achieved by fixing all the parameters except for £1, and AT which are 
varied to minimize the cost function 





(A-15) 



Figure [9] shows the minimization surface profile for the parameters given in Sec. IIVI It is 
important to note that the above cost function leads to a "loose" optimization in the sense 
that it does not strictly optimize the STIRAP into any specific coin operator. Rather, it only 
requires that the coin matrix be unitary. It also turns out that the optimum parameters for 
a unitary C are independent of the choice of phase factors a and (3. Instead these phases can 
be conveniently altered to manipulate the exact form of the operator C while maintaining 
its unitarity. 
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VI. FIGURES 
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FIG. 1: a) An array of quantum dots representing the discrete nodes for a quantum walk on a line, 
b) The electronic structure of a pair of quantum dots QPyiraik Q^Aux- ^^^^ energy 

eigenstates, the overlap between the electron wavefunction inside the dots is negligible. 
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FIG. 2: Labeling of the energy levels as coin states | J,) and | |), excited state |e) and auxiliary 
state \A). STIRAP operations between states | [) < — > \A) and | |) < — > |^) are facilitates via the 
intermediary state |e), using laser pulses with angular frequencies Oj, Of and ^Ia- 
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FIG. 3: Energy level diagram for a single node corresponding to the quantum dot dimensions 
presented in Fig. [1} Ea ~ 173 meV, Ef, ~ 1045 meV and E^ ~ 1912 meV. Energy levels E^ and 
E^ are nearly equidistant from Ea with an energy gap bE^ ~ 15 meV. Similarly, the immediate 
levels above and below Ef. and E^ are separated by bEf, ~ 20 meV and bE^ ~ 30 meV respectively. 
The absorption spectrum is assumed to have a line width ^ 1 meV. Solid lines represent the 
coupling between the desired energy levels via Oj, Vl^ and pulses. Dotted lines demonstrate 
that these frequencies are prevented by the energy structure from activating any spurious coupling 
between other any levels. 
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FIG. 4: The sequence of 2-photon and 3-photon STIRAP to implement a single step of the quantum 
walk, a) The initial state of the walk with the electron confined to state | |, 1); b) A 3-photon 
STIRAP implements the coin rotation C, mixing the states | t)^) ^ind | J,,i); c) A 2-photon 
STIRAP transfers the population from state | to state |j4,i); d) Even barriers are lowered 
and odd barriers are raised in order to regroup the quantum dots. It is now possible for another 
2-photon STIRAP to transfers the population from state \A,i) to state | |,i + 1), completing the 
translation operation t|^; e) Potential barriers are returned to their initial setting and the above 
process repeated. 



16 




Time (ps) 



FIG. 5: (upper panel) The two laser pulses involved in a 2-photon STIRAP process. Both pulses 
have a Gaussian envelope and are applied in the counter intuitive order, i.e. the Stokes pulse S 
responsible for the |e) < — > \A) transition is applied before the pump pulse responsible for the 
I I) < — > |e) transition, (lower panel) The time evolution of dressed states | |) (solid), |e) (dotted) 
and \A) (dashed) due to the application of the 2-photon STIRAP with pulse parameters optimized 
to perform a swap operation. Initially ■0t = 1 and V'A = ipe = 0. 
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FIG. 6: (upper panel) The six laser pulses involved in a 3-photon STIRAP process. All pulses have a 
Gaussian envelope and are applied in the counter intuitive order, (lower panels) The time evolution 
of dressed states | t) (solid), | j) (dashed), |e) (dot-dashed) and \A) (dotted) with pulse parameters 
optimized to implement the coin operators C('^/4, ^2, V^) C{Tr/6, '^/2, ^/2) respectively. Initially 
-i/^l = 1 and ipi = ipA = V'e = 0. 
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FIG. 7: The electron probability distribution after 100 applications of the pulse sequence (solid) 
vs. the corresponding discrete-time quantum random walk distribution after 100 steps (dotted). 
The pulse 3-photon pulse parameters were optimized to perform C(7r/4, '^/2, ^^72) and the electron 
was initially confined to node 1 with probability distribution: = ip^ = (upper panel) and 
^/J^ = 1 and ip^ = (lower panel). 
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FIG. 8: The impact of introducing some noise in the applied pulse parameters, a) Deviation from 
the exact quantum random walk distribution (dashed) due to an induced 2% uncertainty in the 
laser pulse peak energies (solid), b) Deviation from the exact quantum random walk distribution 
(dashed) due to an induced 5% uncertainty in the laser pulse phases (solid), c) Deviation from the 
exact quantum random walk distribution (dashed) due to an induced 2% uncertainty in the laser 
pulse standard deviations (solid), d) Deviation from the exact quantum random walk distribution 
(dashed) due to an induced 0.3% uncertainty in the laser pulse timing (solid). 
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FIG. 9: The optimization surface of the C operator obtained by minimizing the cost function kc 
(Eq. lA-lSp . Referring to Fig. [6l £h is the PI pulse energy peak and AT is the time between the 
S and PI pulse energy peaks. 
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